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Abstract 

In this report we investigate the macroscopic quantum tunnehng of a Bose 
condensate falhng under gravity and scattering on a Gaussian barrier that 
could model a mirror of far-detuned sheet of light. We analyze the effect of 
the inter-atomic interaction and that of a transverse confining potential. We 
show that the quantum tunneling can be quasi-periodic and in this way one 
could generate coherent Bose condensed atomic pulses. In the second part of 
the report, we discuss an effective ID time-dependent non-polynomial nonlin- 
ear Schrodinger equation (NPSE), which describes cigar-shaped condensates. 
NPSE is obtained from the 3D Gross-Pitaevskii equation by using a varia- 
tional approach. We find that NPSE gives much more accurate results than 
all other effective ID equations recently proposed. 

I. INTRODUCTION 

A macroscopic signature of quantum properties of matter is the tunneling of a many- 
particle Bose condensate through a barrier. This subject has been investigated by various 
authors [1-2], in particular in connection to the formation of a Josephson current between 
two wells separated by a potential barrier [3-6]. 



The case we study is that of a faUing condensate that scatters on a potential barrier that 
could model a mirror formed by a far-detuned sheet of hght [7] . We find that the interatomic 
interaction and the geometrical aspects of the system, like the aspect ratio of the cloud or the 
fact the the cloud remains trapped in the transverse directions, have a strong effect on the 
tunneling probability Moreover, we show that in our system macroscopic quantum tunneling 
(MQT) is a quasi-periodic phenomenon and it can be used to generate Bose condensed 
atomic pulses [8]. Finally, from the 3D Gross-Pitaevskii equation and using a variational 
approach, we derive an effective ID wave-equation that describes the axial dynamics of a 
Bose condensate confined in an external potential with cylindrical symmetry. The trapping 
potential is harmonic in the transverse direction and generic in the axial one. Our equation, 
that is a time-dependent non-polynomial nonlinear Schrodinger equation (NPSE), can be 
used to model cigar-shaped condensates. In the limiting cases of weak and strong interaction, 
our approach gives rise to Schrodinger-like equations with different polynomial nonlinearities. 



II. MQT OF A FALLING BOSE CONDENSATE 

The Bose condensate of a dilute gas at zero temperature is well described by the 3D 
Gross-Pitaevskii equation (3D GPE) 



-^V' + VU^)+g{N-imr,t)\' 



V'(r,i), (1) 



where ip{r,t) is the order parameter of the condensate in an external potential Vext{i^), and 
g — Anh^as/m, with the s-wave scattering length. N is the number of condensed atoms. 

The specific system we consider is a condensate falling under gravity and scattering on 
a Gaussian potential barrier. At t < the condensate is confined by a harmonic trap and 
the external potential reads 

TTl ^ 

Vextip, z) = -^{^y + ^l{z- zof) + mgz + Ue'^ , (2) 

where p = (x^ + y^Y^"^ . We use harmonic oscillator units: ujh = {'-^p'-^zY^^ = 2% x 100 Hz. 
^^Na atoms with an — {h/{'>TUVH)y^'^ — 27 //m and = 3 nm. 



We set zq — ISoij so that the condensate is initially far from the Gaussian potential 
barrier. The 3D GPE is numerically solved by using a predictor-corrector splitting method 
[8] . We calculate the ground-state wave- function in the external potential by means of the 
splitting method with imaginary time. Then, we switch off the harmonic potential and 
use the previous wave-function as initial condition. The total energy per particle of the 
condensate is about 180 hurn with N ranging from 1 to 10^. 

For large values of the energy barrier there is pure bouncing with interference. This 
effect have been experimentally observed by Bongs et al [7] . By reducing the energy barrier, 
in addition to bouncing there is MQT. We have found that, if initially the cloud is spherical 
under free fall, the tunneling fraction grows with the number of particles, i.e. the chemical 
potential. Moreover, as shown by Fig. 2, MQT is quasi-periodic phenomenon due to the 
quasi-periodic bouncing of a reflected part of the condensate, that is also expanding. Such 
a mechanism could be used to generate Bose condensed atomic pulses. 

Fig. 3 shows that the tunneling fraction is reduced by the interatomic interaction if 
the initially spherical falling condensate is under a strong enough transverse conflnement. 
Instead, for a cigar-shaped quasi- ID Bose condensate the interatomic interaction enhances 
the tunneling rate, in agreement with the theoretical predictions of Ref. [3]. 

III. EFFECTIVE ID EQUATION FOR CIGAR-SHAPED CONDENSATES 

As previously stated, the 3D Gross-Pitaevskii equation (3D GPE) is accurate to describe 
a Bose condensate of dilute bosons near zero temperature, where thermal excitations can 
be neglected. An interesting problem is to find a reliable mapping from the 3D GPE to an 
effective ID equation that describes cigar-shaped Bose condensates. This problem is non 
trivial due to the nonlinearity of the GPE. In this section we illustrate a new variational 
approach which gives very accurate results. 

First we observe that the 3D GPE can be obtained by minimizing the following action 
functional 
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Note that we now use TV instead of — 1. We analyze a trapping potential that is generic 
in the axial direction and harmonic in the radial one: 



V,,t{v) = -mujl{x^ + y^) + V{z) (4) 
We minimize the action functional S by choosing the following trial wavefunction [9] : 

V^(r, t) = (i){x, y, t; r]{z, t)) f{z, t) , (5) 
where both and / are normalized and 4> is represented by a Gaussian [10]: 

4>{x,y,t-r]{z,t))^^^^^j^^^. (6) 

The variational functions r){z,t) and f{z,t) will be determined by minimizing the action 
functional after integration in the {x, y) plane. 

From the two Euler-Lagrange equations, one eventually obtains 
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This equation is a time-dependent non-polynomial nonlinear Schrodinger equation (NPSE) 
[10]. Note that from NPSE in certain limiting cases one recovers familiar results. In the 
weakly-interacting limit agA^j/p << 1 the previous equation reduces to 
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where the constant hu!± has been omitted because it does not affect the dynamics. This 
equation is a ID Gross-Pitaevskii equation (ID GPE), whose nonhnear coefficient is g' — 

g/{2Tia\_). This ansatz has been used by various authors, for example in Ref. [11]. 
In the strongly-interacting limit asN\f\'^ » 1 the NPSE becomes 



/• (9) 
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In this limit, and in the stationary case, the kinetic term can be neglected (Thomas- Fermi 
approximation) and one finds 

where 11' is the chemical potential, fixed by the normalization condition. This ID Thomas- 
Fermi density profile is quadratic in the term fi' — V{z). The same quadratic depencence 
is obtained starting from the Thomas-Fermi approximation of the 3D stationary GPE, i.e. 
neglecting the spatial derivatives in Eq. (1), and then integrating along x and y variables. 
In this way one finds a formula that difi^ers from the previous one only for the numerical 
factor which is 1/4 instead of 2/9. 

To test the accuracy of the NPSE, we take V{z) — {l/2)muj1z^ , and compare the results 
of 4 different procedures: 1) 3D GPE: numerical solution of 3D GPE; 2) ID GPE: numerical 
solution of ID GPE. 3) CGPE: numerical solution of ID GPE, where the nonlinear term 
is found by imposing that the ID wave- function has the same chemical potential of the 3D 
one [12]; 4) NPSE: numerical solution of NPSE. 

As shown in Fig. 5, CGPE gives better results than ID GPE for large values of scattering 
length but in any case NPSE is superior. Our calculations (for details see [10]) show that the 
ground-state and also the dynamics of the condensate are accurately described by NPSE. 



CONCLUSIONS 

We have studied the dynamics of Bose condensates falling under gravity and scattering on 
a Gaussian barrier that models a mirror of light. Apart the pure bouncing with interference, 
that can be seen for very large values of the energy barrier, we have investigated quantum 
tunneling. In our system quantum tunneling is a quasi-periodic phenomenon: it can be 
used to generate Bose-Einstein condensed atomic pulses. In addition, we have found that 
geometrical aspects are important in tunneling: ID results do not always reflect the behavior 
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of 3D systems. Finally, we have verified that the ID non-polynomial nonlinear Schrodinger 
equation we have obtained by using a Gaussian variational ansatz from the 3D Gross- 
Pitaevskii action functional is quite reliable in describing cigar-shaped condensates. The 
agreement with the results of the 3D Gross-Pitaevskii equation is very good for both ground- 
state and dynamics of the condensate. This equation will be useful for detailed numerical 
analysis of the dynamics of cigar-shaped condensates, particularly when the local density 
may undergo sudden and large variations. In conclusion, we observe that Pulsed MQT of 
matter waves can be generated also with other mechanisms. For example, we are currently 
studying the stimulated MQT of matter waves produced by a Bose condensate trapped in 
a potential with a periodically modulated energy barrier. 
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Fig. 1: Density profile = 0, z)^ of the free falling Bose condensate at t — 0, t — 2 

1 I 2 

and t = 4. Gaussian energy barrier Ue^ ''^ , with U = 200 and u = 1. Initial chemical 
potential: ii = 184.23. N = 10^ condensed atoms. Lenght in units an = 27 /xm, time in 
units ujj^ = 1.6 ms and energy in units huju- 
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Fig. 2: Tunneling fraction as a function of time t. Comparison between radially-free 
{(jp — 0) and under radial confinament {cOp = 0.5) falling condensates. N — 10^ condensed 

2/2 

atoms and initial position zq = 15. Gaussian energy barrier Ue^ ''^ , with a — 1. Units as 
in Fig 1. 
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Fig. 3: Tunneling fraction as a function of the initial chemical potential /j, for different 
frequencies cUp of radial confinament. The points in each line correspond to a sequence of 
numbers N of atoms. From left to right: N = 1, 10^, 10^, 10^, 5 x 10^, 10^. Gaussian energy 
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barrier Ue^ ''^ , with U = 200 and a = 1. Units as in Fig 1. 
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Fig. 4: Density profile \ijj{p = 0, at t = for various values of the number N of 

2 / 2 

condensed atoms. Gaussian energy barrier Ue^ , with U — 200 and a — 1. Radially- 
confined {ujp = 1) falling condensate with cigar-shaped {cUp/cUz = 10) initial condition. Pp is 
the tunneling fraction. Units as in Fig 1. 
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Fig. 5: Normalized density profile p{z) = along the axial direction z of the 

Bose condensate in harmonic potential. Number of Bosons: N — 10^ and trap anisotropy: 
'^lI'jJz = 10. Four different procedures: 3D GPE (solid line), ID GPE (dashed line), 
CGPE (long-dashed line) and NPSE (dotted line). From top to bottom: a^/a^ = 10~^, 
Os/oz = 10~^, tts/az — 10~^. Length z in units of az — Jh,/{mujz)- 
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